Charge Densities and Charge Noise in Mesoscopic Conductors 
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We introduce a hierarchy of density of states to characterize the charge distribution in a meso- 
scopic conductor. At the bottom of this hierarchy are the partial density of states which represent 
the contribution to the local density of states if both the incident and the out-going scattering 
channel is prescribed. The partial density of states play a prominent role in measurements with a 
scanning tunneling microscope on multiprobe conductors in the presence of current flow. The partial 
density of states determine the degree of dephasing generated by a weakly coupled voltage probe. 
In addition the partial density of states determine the frequency-dependent response of mesoscopic 
conductors in the presence of slowly oscillating voltages applied to the contacts of the sample. The 
partial density of states permit the formulation of a Friedel sum rule which can be applied locally. 
We introduce the off-diagonal elements of the partial density of states matrix to describe charge 
fluctuation processes. This generalization leads to a local Wigner-Smith life-time matrix. 



The characterization of the charge distribution and its 
fluctuations are a central aspect of electrical conduction 
theory. In this work we are concerned with mesoscopic 
conductors, structures which are so small and cooled to 
such low temperatures that the wave nature of electron 
motion becomes apparent. Three recent reviews Bijovidc 
an entry to the literature of this very active fielcaia. Of 
interest is the charge distribution in such a conductor and 
its fluctuations under equilibrium and non-equilibrium 
conditions. The charge distribution is an intersting quan- 
tity in itself. However, in addition, it is important for the 
solution of a number of problems. We mention here only 
a few: Already, in the equilibrium state, the charge dis- 
tribution needs to be known to determine the effective 
equilibrium one-electron potential. It is this potential 
which determines the scattering properties of the con- 
ductor and thus the conductance coefficients. If volt- 
ages are applied which drive the conductor away from 
the equilibrium state, typically, the charge distribution 
also undergoes a bias dependent modification. Thus the 
charge distribution is needed to find the non-linear I- 
V-characteristic of an electrical conductor. Nanoscopic 
contacts to a sample measure local density of states and 
thus a charge distribution. Knowledge of the charge dis- 
tribution is essential for the discussion of the dynamic 
conductance like the ac-conductance or photon assisted 
transport. The charge distribution also determines the 
capacitance coefficients of various spatial regions within 
the mesoscopic structure and also vis-a-vis nearby gates. 



I. LOCAL PARTIAL DENSITY OF STATES 

It is the purpose of this_work to present a discussion 
of the charge distributioroEl and the charge noiseH based 
on the scattering .approach. Similar to the discussion of 
linear transportQcl, we view the sample as a target which 
transmits or reflects electrical carriers into the contacts 
connected to the sample. A particular examplc£3 of a 



structure of interest here is shown in Fig. 1. For purely 
elastic scattering the sample is described by a scattering 
matrix s. We label the different contacts to the sample 
by a — 1,2, .... Each contact is treated as perfect wave 
guide which permits the definition of incoming and out- 
going quantum channels. Each contact has at a given 
energy a set of quantum channels with quantum num- 
bers labeled m. The scattering matrix element s am ^p n 
relates the out going current amplitude in contact a in 
channel m to the incoming current amplitude in contact 
(3 in channel n. For the following discussion it is useful to 
introduce the submatrice s Q( g which combines all scatter- 
ing matrix elements with the same contact indices. This 
is a matrix of dimension M a xMp for a conductor with 
M a open scattering channels in contact a and Mp open 
channels in contact j3. We will now show that the equi- 
librium and non-equilibrium charge distribution can be 
described with the help of a hierarchy of density of states. 
On top of this pyramid is the local density of states. At 
the bottom nf this pyramid are the local partial density 
of stateMW 
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which give the contribution to the local density of states 
at point r of carriers incident in contact (3 and leaving 
the conductor through contact a. In Eq. (Q) the trace 
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FIG. 1. Mesoscopic conductor with a tunneling contact 
with coupling energy t &t x. After Ref. [j^. 
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represents a summation over channel indices. In the par- 
tial density of states we specify both the incident con- 
tact (a pre-selection) and the exiting contact (a post- 
selection). The local partial density of states is obtained 
by investigating the change of the scattering matrix due 
to a small additional potential perturbation SU(r) at 
point r. 

On the next higher level in the pyramid of density of 
states are the injectivities and emissivities. The injec- 
tivity of contact j3 gives the contribution to the density 
of states of carriers entering the conductor through (3 ir- 
respective of the contact through which they leave the 
sample. The injectivity i/(r, 0) into point r of contact 
(3 is obtained by summing all partial density of states 
via, r, f3) over the first contact index, 
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Using the unitarity of the scattering matrix we can also 
write for the injectivity, 
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The emissivity of point r into contact a is that portion 
of the local density of states which consists of carriers 
which leave the sample through contact a irrespective 
of the contact through which they entered the sample. 
It is obtained by summing the partial density of states 
over the second contact index (3. Taking into account the 
unitarity of the scattering matrix s we find, 



v(a,r) 



1 

2tt7 



13 



Ss a 



/3 



"i3 



eSU(r) 



(4) 



At the top of this pyramid of density of states is the local 
density of states. For the local density of states we im- 
pose no restriction as to where the carriers come from nor 
as to where they will go in the future. The local density 
of states is thus the sum over all local partial density of 
states. In particular it is also the sum of all injectivities 
or the sum of all emissivities, 



and is thus given by 
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We will later show the usefulness of these density of 
states. Here we only briefly mention the magnetic field 
symmetry. Since in a uniform field the scattering ma- 
trix has the reciprocity symmetry s a ^(_B) = s a p(—B), 
we see that the partial density of states have the sym- 
metry vb{ol, r, 0) = v_b(P> r , <*). Similarly, the injectivi- 
ties and emissivities are related by reciprocity ^s(r, a) — 



i/_B(a, r) whereas the local density of states is an even 
function of magnetic field ^s(r) = u-b{y). Fig. 2 shows 
the local density of states, the injectivity and the ratio of 
the injectivity to the local density of states for the sim- 
ple case of a delta function barrier in an otherwise per- 
fect ballistic one-channel wire connecting contacts 1 and 
2. For additional simple examples we refer the reader 
to Refs. p~0| - |T^ . Partial density of states in disordered 
structures are evaluated numerically in Ref. \LSi 
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FIG. 2. The local density of states v(x) (solid line), the 
injectivity v{x, 1) of contact 1 (dotted line) (both in units 
of 1/hv), and the ratio of these two densities v(x,l)/v(x) 
(dashed line) of a one dimensional wire with a 8 barrier at 
x — which a strength that leads to a transmission proba- 
bility of T = 0.8. The coordinate x is in units of the Fermi 
wavelength A. After Ref. 



II. GLOBAL PARTIAL DENSITY OF STATES 

Often we are interested not in the density of states at 
a given point but in a certain volume O. Of course the 
partial density of states in f2 are obtained simply by inte- 
grating the local partial density of states over the volume 
of interest. Of particular interest are the density of states 
of the entire sample. In this case the surface of the vol- 
ume fi will intersect the leads connected to the sample. 
Suppose now that these intersections coincide with our 
specification of the scattering matrix. The phases which 
occur in the scattering matrix s a p are the phases which 
are accumulated through scattering between the inter- 
section S a of CI in lead a and the intersection of Sp in 
lead (3. Under these circumstances a simplification of the 
above expressions is possible, since the functional deriva- 
tive with respect to the local potential and the integration 
over f2 can be replaced by an energy derivative, 
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For the global density of states (density of states of the 
entire sample) we find the partial density of states, 
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from which we obtain the injectance of contact a 
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the emittance into contact a. 
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and the global density of states 
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The scattering matrix is unitary and all its eigenvalues 
A, = exp(iQ) are thus on the unit circle. In terms of 
these eigenvalues we have 
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where det(s) is the determinant of the scattering matrix. 

Eq.j— i(M) is a familiar expression of the density of 
stateslij El We recall it here to emphasize that it is 
an approximate expression for the density of states, o£ a 
mesoscopic conductor. In nuclear scattering theorjJii'Eil 
it is turned into an exact expression by considering not 
the total density of states like we have done but the ex- 
cess density of states in comparison to a free problem and 
by taking the volume of interest to infinity (The problem 
with and without a scattering object are compared). In 
contrast for our type of scattering problems (for instance 
for a multiprobe mesoscopic conductor) there is no rea- 
sonably defined free problem and similarly it would be 
meaningless to consider the volume of interest to become 
infinitely large. Therefore in mesoscopic physics theMx*. 
act expression for the density of states of the sample iafrlLil 
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whereas the more familiar expression Eq. (12) is an ap- 
proximation. 

Except in very special cases the two expressions for 
the local density need not to be distinguished. However 
the difference is important when WKB does not hold. 
For instance using Eq. (p"2|), the density of states in- 
side a rectangular tunneling barrier of height Uq diverges 
like i? -1 / 2 as the energy kinetic energy E of the incident 
particles tends to zero whereas if one uses Eq. ({l3|) one 
finds (the much more reasonable result) that the den- 
sity of states is finite and vanishes like E 1 / 2 for small 
energiesO. The distinction between energy derivatives 
and potential derivatives is fam Uiai . fr am the discussion 
of characteristic tunneling timeal3'EJ'llac2l. For a discus- 
sion of the connection between the local partial density 
of states and characteristic tunneling times we refer the 
reader to Ref. Ell 



III. ABSORPTION AND EMISSION OF 
PARTICLES: INACTIVITIES AND 
EMISSIVITIES 



We next present a simple example which illustrates 
the usefulness of the densities of states introduced above. 
Suppose that the potential has in same small volume ele- 
ment fl a small imaginary partE3 cj. Thus the potential 
is £7(r) out-side but is U(r) — iUT inside f2. To pro- 
ceed we take into account that the scattering matrix is 
a functional of the potential and expand the scattering 
matrix to first order in T. This gives for the transmis- 
sion and reflection probabilities in the presence of a small 
absorptior£3 



T a0 (T) = Ta/3 ^1 - dr T u(a,T,0) 



(14) 



where T a /3 is the transmission probability of the scatter- 
ing problem without absorption if a and (3 are different. 
If a and (3 are equal, we find a reflection probability 



-Raa(r) = i?aa ^1 — J dr T v(a, i", a) 
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Thus in the presence of absorption the reflection and 
transmission probabilities are " diminished" according to 
the partial density of states. Consider now the carriers 
streaming from contact j3 into the sample in a narrow en- 
ergy interval dE. The incident current Ii n must be equal 
to the sum of the transmitted current It, the reflected 
current In and the absorbed current It, 



It+Ir + It- 



(16) 



Taking into account that the total incident current from 
contact a is Ip = {e/h)MpdE, and taking into account 
that J2aj£t3 T a f3 + Rpp = Mp we find that the total 
abosrbed flux is given by the sum 
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Thus the absorbed flux is determined by the injectivity 
of contact (3 into the volume of interest. 

Another way of determining the absorbed flux proceeds 
as follows. The absorbed flux is proportional to the in- 
tegrated density of particles in the region of absorption. 
The density of particles can be found from the scattering 
state ipp m (r) which describes carriers incident from con- 
tact f3 in channel m. The wave function is normalized 
such that it has unit amplitude in the incident channel. 
We assume that all incident channels m. contact (3 are 
filled equally. The absorbed flux is thusE3 
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Note that here the density of states 1 / hvp m of the asymp- 
totic scattering region appears. vp m is the velocity of 
carriers in the incident channel. Thus we have found a 
wave function representation for the injectivity. Compar- 
ing Eq. (pi) and Eq. (Ill) gives 
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The total local density of states v(y) at point r is obtained 
by considering carriers incident from both contacts. 

There is now an interesting additional problem to be 
addressed. Instead of a potential which acts as a carrier 
sink (as an absorber) we can ask about a potential which 
acts as a carrier source. Obviously, all we have to do to 
turn our potential into a carrier source is to change the 
sign of the imaginary part of the potential. As a con- 
sequence of a carrier source in the volume f2 we should 
observe a particle current into all contacts of our sample. 
To find these currents we consider a small energy interval 
at the energy of interest. The total current injected into 
the sample in Q is 



hn(y) = / drTi/(r). 
Jn 



(20) 



Taking into account that the incident current is normal- 
ized to 1 the current I ou t (/?) in contact /? due to a carrier 
source in Q is given by 

W") = {e/h)dE[M a - R aa (T) - £ T aB (T)] (21) 

due to the modification of both the transmission and re- 
flection coefficients. Using Eqs. ( |l4j ) and jl5| ) (with T 
replaced by -r) gives 

lout («) = - / dr T V v(a, r,0) = - [drT v{a, r). 
Jn g Jn 

(22) 

The current in contact a is determined by the emissivity 
v(a, r) of the region Q into contact a. 

We conclude this section with a remark of caution. The 
partial densities are not necessarily positive as we would 
expect from a usual density of states. (The simple ex- 
ample of an (asymmetric) resonant double barrier shows 
that one of the diagonal elcmcpts v{a, r, a) has a range 
of energies where it is negativeca Note that this implies 
that a small absorption inside a sample can increase a re- 
flection probability) . The fact that the partial density of 
states are not necessarily positive is a consequence of over 
specification: both the incident and the exiting channel 
are prescribed. In contrast, the injectivities and emissiv- 
ities are, however, always positive. The proof is given by 
Eq. @. 



IV. GENERALIZED BARDEEN FORMULAE 



It is well known that with a scanning tunneling mi- 
croscopje (STM) we can measure the local density of 
statesEj. STM measurements are typically performed in 
a two terminal geometry, in which the tip of the mi- 
croscope represents one contact and the sample provides 
another contacto. Here we are interested in the trans- 
mission probability from an STM tip into the contact of 
a sample with two or more contacts as shown in Fig. |]. 
Geometries which—use more than one tunneling tip are 
also of interestE£lE3. Thus we deal with a multiterminal 
transmission problemEij. If we denote the contacts of the 
sample by a Greek letter a = 1,2,.. and use tip to la- 
bel the contact of the STM tip, we are interested in the 
tunneling probabilities T at i P from the tip into contact a 
of the sample. In this case the STM tip acts as carrier 
source. Similarly we ask about the transmission prob- 
ability Ttip a from a sample contact to the tip. In this 
case the STM tip acts as a carrier sink. Earlier work 
has addressed thisHpeoblem either with the help of elec- 
tron .wave dividersE3'E3, or by applying the Fermi golden 
Rulea. Recently, Gramespacher and the authorEHl have 
returned to this problem and have derived expressions 
for these transmission probabilities fxom the scattering 
matrix of the full multiprobe problemO (sample plus tip) . 
For a tunneling contact with a density of states vtip which 
couples locally at the point x with a coupling energy \t\ 
these authors found 



Tti P ,a = ^ 2 v t ip\t\ 2 v{x,a) , 



T a up = ^TT 2 v(a,x)\t\ 2 v tip . 



(23) 



(24) 



In a multiterminal sample the transmission probability 
from a contact a to the STM tip is given by the injec- 
tivity of contact a into the point x and the transmission 
probability from the tip to the contact a is given by the 
emissivity of the point x into contact a. Eqs. (23) and 
( p4| ) when multiplied by the unit of conductance e 2 /h are 
generalized Bardeen conductances for multiprobe con- 
ductors. Since the local density of states of the tip is 
an even function of magnetic field and since the injectiv- 
ity and emissivity are related by reciprocity we also have 
the reciprocity relation Tu p>a (B) — T a . t ip(—B). 

The presence of the tip also affects transmission and 
reflection at the massive contacts of the sample. To first 
order in the coupling energy \t\ 2 these probabilities are 
given by 



rptip 

1 af3 
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(25) 



Here the upper index tip indicates that this transmission 
probability is calculated in the presence of the tip. T a p 
without an upper index is the transmission probability 
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in the absence of the tip. The correction to the trans- 
mission probabilities a ^ f3 and reflection probabilities 
a = (3 is determined by the partial densities of states, 
the coupling energy and the density of states in the tip. 
Note that if these probabilities are placed in a matrix 
then each row and each column of this matrix adds up 
to the number of quantum channels in the contacts. 



V. VOLTAGE PROBE AND INELASTIC 
SCATTERING 

Consider a two probe conductor much smaller than any 
inelastic or phase breaking length. The conductance of 
such a conductor can then be said to be coherent and 
its conductance is at zero temperature given by the Lan- 
dauer formula G = (e 2 /h)T, where T is the probability 
for transmission form one contact to the other. How is 
this result affected by events which break the phase or 
by events which cause inelastic scattering? To investigate 
this question Ref. |3FJ proposes to use an additional (third) 
contact to the sample. The third probe acts as a voltage 
probe which has its potential adjusted in such a way that 
there is no net current flowing into this additional probe, 
Is = 0. The current at the third probe is set to zero by 
floating the voltage /13 = eV$ at this contact to a value 
for which ^3 vanishes. The third probe acts, therefore, 
like a voltage probe. Even though the total current at the 
voltage probe vanishes individual carriers can enter this 
probe if they are at the same time replaced by carriers 
emanating from the probeEH Entering and leaving a con- 
tact are irreversible processes, since there is no definite 
phase relationship between a carrier that enters the con- 
tact and a carrier that leaves a contact. In a three probe 
conductor, the relationship between currents and volt- 
ages is given by I a = J2p G^Vp where the G a n are the 
conductance coefficients. Using the condition ^3 = to 
find the potential V3 and eliminating this potential in the 
equation for I2 or 1\ gives for the two probe conductance 
in presence of the voltage probe 



G = -{G 
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(26) 



For a very weakly coupled voltage probe (see Fig. [I]) 
we can use Eqs. (|23| - |25|). Taking into account that 
G a p = -(e 2 //i)T*i for a ^ (3 we find 



G=-(T-47r 2 |i|>(2,x,I) 



v(2,x)v(x, I) 
v(x) 



])• (27) 



Here v(x) is the local density of states at the location of 
the point at which the voltage probe couples_to the con- 
ductor. Eq. (27) has a simple interpretatiortJ. The first 
term T is the transmission probability of the conductor 
in the absence of the voltage probe. The first term inside 
the brackets proportional to the local partial density of 
states gives the reduction of coherent transmission due 



to the presence of the voltage probe. The second term in 
the brackets is the incoherent contribution to transport 
due to inelastic scattering induced by the voltage probe. 
It is proportional to the injectivity of contact 1 at point 
x. A fraction v(2 1 x)/v{x) of the carriers which reach 
this point, proportional to its emissivity, are scattered 
forward and, therefore, contribute to transport. Notice 
the different signs of these two contributions. The effect 
of inelastic scattering (or dephasing) can either enhance 
transport or diminish transport, depending on whether 
the reduction of coherent transmission (first term) or the 
increase due to incoherent transmission (second term) 
dominates. 

Instead of a voltage probe, we can also use an opti- 
cal potential to simulate inelastic scattering or dephas- 
ing. However, in order to preserve current, we must use 
both an absorbing optical potential (to take carriers out) 
and an emitting optical potential (to reinsert carriers). 
The absorbed and re^emitted current must again exactly 
balance each otherEll. From Eq. (|I^) it is seen that 
the coherent current is again diminished by Yv(2,x, 1), 
i. e. by the partial density of states at point x. The 
total absorbed current is proportional to Yv(x,\\ the 
injectance of contact 1 into this point. As shown in 
section 3 a carrier emitting optical potential at x gen- 
erates a current — Tu(l,x) in contact I and generates a 
current —Tv(2, x) in contact 2. It produces thus a to- 
tal current — Tv[x). In order that the generated and 
the absorbed current are equal we have to normalize 
the emitting optical potentialcil such that it generates 
a total current proportional to Tv(x, I). The current 
produced by an optical potential normalized in such a 
way is thus -pTz/(2, x)v(x, \)/v{x). The sum of the two 
contributionstJ, the absorbed current and the re-emitted 
current gives an overall transmission (or conductance) 
which is given by Eq. (27) with 47r|i| 2 replaced by T. 



Thus the weakly coupled voltage probe (which has cur- 
rent conservation built in) and a discussion based on 
optical potentials coupled with a current conserving re- 
insertion of carriers are equivalentEj. There are discus- 
sions in the literature which invoke optical potentials but 
do not re- insert carriers. Obviously, such discussions vi- 
olate current conservation. Even if current is conserved, 
the partitioning of the reinserted carriers has to be done 
as described above. Otherwise such fundamental proper- 
ties like the Onsager symmetry are violatecu3. 

We would like to point out that there is in general 
no reason that tunnel contacts and absorption and re- 
emission with the help of imaginary potentials lead to 
identical results. The equivalence rests on a particular 
description of the voltage probe. There are many differ- 
ent ways of coupling a tunnel contact and our description 
of the voltage probe given here is not unique. The claim 
can only be that for sufficiently weak optical absorption 
and re-insertion of carriers there exits one voltage probe 
model which gives the same answer. 
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VI. THE FRIEDEL SUM RULE 

The connection between the phase shift of scattering 
elements and the charge is known as Friedel sum rule. 
It is most often stated as a global relation between the 
charge in an energy interval dE and the derivatives of the 
scattering phases with respect to energy as 



dE 2vr dE 



log[det(s)]. 



(28) 



From the discussion at the end of Section II it is evident 
that this relation is not exact since in the applications 
envisioned here the scattering matrix and the phases are 
calculated for a finite volume fl. Thus whenever WKB 
does not hold Eq. (|2lj| ) will not hold. The exact relation 
follows from Eq. (|19|) and Eqs. (||) and (||) and is given 
by 
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This relation is exact locally and by integration can be ap- 
plied to any volume of interest. Application of Eq. (|2S| ) 
to situations in which WKB does not hold has already 
given rise to statements in the literature that the Etiedel 
sum rule does not apply in mesoscopic conductorsE3. 



VII. THE CHARGE RESPONSE TO AN 
INTERNAL POTENTIAL 

We next discuss problems for which the self-consistent 
adjustment of the potential is important. With this goal 
in mind we first investigate the variation of the local car- 
rier density Sn(r) in response to a local potential vari- 
ation 6U(r') at a point r'. This variation is given by 
the Lindhard function, Sn(r) = — J dr' H(r,r')eSU(r'). 
This function plays an important role if we want to de- 
termine the carrier distribution self-consistently. Since 
the local charge density at equilibrium is n(r) = 
J dEv(r)f(E), it follows immediately that the Lindhard 
function is given by 



n(r, r') 



dE 



Sv(t) 
e5U{r>) 



f(E). 



(30) 



This response now involves second order functional 
derivatives of the scattering matrix. It can be simpli- 
fied, if we are interested in potential variations not on 
the microscopic scale but only over distances large to a 
Fermi wave length Xp. In that case we can integrate 
the Lindhard function over its second argument over a 
volume that is large compared to a Fermi wave length 



and use the WKB formula which replaces the functional 
derivative with respect to U(r') with an energy deriva- 
tive. This gives 

n(r,r') = 5(r-r') J dE^-f(E). (31) 

Partial integration with respect to the energy E gives 

B(r', r) = 5(t> -r)J dEv{v) (-^) • (32) 

Over sufficiently large distances the Lindhard function is 
determined by the local density of states. 



VIII. THE LOCAL CHARGE RESPONSE OF A 
MESOSCOPIC CONDUCTOR 

Consider a small deviation Sfx a = fx a — fj. of the elec- 
trochemical potential fj, a of contact a away from its 
equilibrium value fi. At fixed internal potential this 
increases the carrier density inside the conductor by 
dn(r) = v(r,a)5fj, a . This injected charge generates 
due to the Coulomb interaction an electrostatic response 
SU(r) which via the Lindhard function will also give a 
contribution to the charge density. Thus the local carrier 
density isU 



6n(r) 



>(r,a)5/x 



J dr' 3 Il{r',r)5U(r') 



(33) 



Suppose that all potentials Sfi a applied to the conductor 
are rised by the same amount S/i. Since the state of a 
conductor can only be changed by voltage differences it 
follows that the electrostatic potential increases also by 
SU(r') = S/i and that Sn(r) ff 0. As a consequence of 
this gauge invariance we havea 



v(r, a) = v{t) = 



J rfr' 3 n(r',r). 



(34) 



The spatial integral over one of the arguments of the 
Lindhard function is equal to the local density of states. 
Note that Eq. ( |30| ) is a symmetric function of its spatial 
arguments. 

To bring the determination of the local charge response 
to an end, we must now insert Eq. ( |33"|) into the Pois- 
son equation and solve for the potential SU(r'). Sup- 
pose that a test charge at point r' generates at point 
r a potential g(v, r'). Here g(r,r') is the Green's func- 
tion of the Poisson equation with a non-local screening 
kernel determined by the Lindhard function mentioned 
above. The injected charge ev(r, a)Sfj, a generates the. 
potential M Q (r)<5/x a where the characteristic potentialtj 
u a (r) is given by 



J dr' 3 5 (r,r>(r',a). 



(35) 



G 



The injectivity is the source term of the characteristic 
potential u a (r). With the help of the characteristic po- 
tentials the internal potential can be expressed in linear 
response to the applied electrochemical potentials. Using 
this, the local carrier distribution is found to be 

5n(r) = ^y ot (r,a)<5^a (36) 

a 

with a screened density of states! 

z/°'(r,a) = v{r,a) -J dr' 3 dr" 3 n(r, r').g(r', r'>(r", a). 

(37) 

Eq. ( |36| ) provides a gauge invariant description of the lo- 
cal carrier density in response to small variations of the 
potentials applied to the contacts of the sample. The 
total variation of the charge is zero, J n drSn(r) = 0, for 
a volume f2 that is sufficiently large (encloses the entire 
sample). 



IX. DYNAMIC CONDUCTANCE 

In this section we discuss as an additional application 
of partial density of states briefly the ac- conductance of 
mesoscopic systems. We consider a conductor with an 
arbitrary number of contacts labeled by a Greek index 
a = 1,2,3.... The problem is to find the relationship 
between the currents I a (u>) at frequnecy u) measured at 
the contacts of the sample in response to a sinusoidal 
voltage with amplitude Vp(u>) applied to contact (3. The 
relationship between currents and voltages is given by 
a dynamical conductance matrixcra G a p{uj) such that 
I a (u)) = G a {3(tjj)V{3(uo). All electric fields are local- 
ized in space. The overall charge on the conductor is 
conserved. Consequently, current is also conserved and 
the currents depend only on voltage differences. Current 
conservation implies ^2 a G a p — for each (3. In order 
that only voltage differences matter, the dynamical con- 
ductance matrix has to obey G a p — for each a. 
We are interested here in the low frequency behavior of 
the conductance and therefore we can expand the con- 
ductance in powers of the frequency^, 

G a p(w) = G% - nuE aP + K aP uj 2 + 0(lu 3 ). (38) 

Here G° a is the dc-conductance matrix. E a p is called the 
emittance matrix and governs the displacement currents. 
Kap gives the response to second order in the frequency. 
All matrices G^p, E a p and K a p are real. 

We focus here on the emittance matrix E a p . The con- 
servation of the total charge can only be achieved by 
considering the long-range Coulomb interaction. To first 
order in frequency it is precisely the charge distribution 
given by Eq. ([36]) which counts. Ref. [| finds for the 
emittance matrix 



E af3 = e 2 [J drv{a, r, (3) - J dr'dTu{a, r')g(r', r>(r, (3)] 

(39) 

Here the first term proportional to the integrated par- 
tial density of states is the ac-response at low frequencies 
which wc would have in the absence of interactions. The 
second term has the following simple interpretation: an 
ac-voltage applied to contact (3 would (in the absence of 
interactions) lead to a charge built up at point r given 
by the injectivity of contact f3. Due to interaction, this 
charge generates at point r' a variation in the local po- 
tential which then induces a current in contact a propor- 
tional to the emissivity of this point into contact a. Note, 
that if screening is local (over a length scale of a Thomas 
Fermi wave length) we have <?(r',r) = 5(r' — rW~ 1 (r). 
In this limit the close connection between Eq. (39) and 
Eq. (|24|)r-i¥ then obvious. We have to refer the reader to 
a reviewO for a more extended discussion and for refer- 
ences to specific applications of Eq. ( |39"| ) . 

A closely related development which could be discussed 
here is a theory of quantum pumping in small systems. In 
quantum pumping one is interested in the current gen- 
erated as two parameters (like gate voltages, magnetic 
flux) which modulate the system are varyed sinusoidally 
but out of phase— J3rouwert3, Avron et al.E3, and Polian- 
ski and BrouwcrL3 develop a theory which is based on the 
slow modulation of the partial densities of states. 



X. PARTIAL DENSITY OF STATES MATRIX 

Thus far we have been interested in the charge distribu- 
tion inside a mesoscopic conductor which is established 
in response to changes in external parameters, such as 
a voltage applied to the contact of the sample. Similar 
to the dynamic fluctuations of the current away from its 
average value, we can also ask about the fluctuations of 
the charge away from its average value. Spontaneous dy- 
namic fluctuations of current and charge can have many 
sources in an electrical conductor. Here we are con- 
cerned with contributions which are fundamental in the 
sense that they can not be avoided. Nyquist or thermal 
noise can not be avoided at any elevated temperature. 
A second source, shot noise results from the granular- 
ity of charge (the charge quantization) and the quantum 
mechanical uncertainty which arises whenever there are 
different possible outcomes of a scattering processu (such 
as the possibility of a particle to be either reflected or 
transmitted). 

To describepthe charge fluctuations due to carriers in 
the conductorQ, we consider the Fermi-field 



*(r,t)=X; / dE 



tp am (r,E) 

(hv am (E)y/- 



;a am (E)exp(-iEt/h), 
(40) 
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which annihilates an electron at point r and time t. The 
Fermi field Eq. ([lo]) is built up from all scattering states 
V'amtr, E) which have unit incident amplitude in contact 
a in channel m. The operator a am (E) annihilates an 
incident carrier in reservoir a in channel m. v am is the 
velocity in the incident channel m in reservoir a. The 
local carrier density at point r and time t is determined 
by n(r, t) = ^'(r, ^^(r, t). We will investigate the den- 
sity operator in the frequency domain, n(r, to). Using the 
Fermi-field we find, 



n(r, uj) 



E 

am{3n 



dE 



1 



1 



(hv am (E)y/i (hv 0n (E + htu)y/2 



the average charge density. Inserting Eq. into the 

Poisson equation with the non-local screening kernel de- 
termines a potential operatora 

&ap(T,w) = J dr , g( r y)al(E)v a f l (i')kf l (E + hu). (46) 

In random phase approximations the potential fluctua- 
tions give rise to induced charge fluctuations. Thus the 
true low frequency charge fluctuations in a conductor are 
determined by the contribution of the bare charge fluc- 
tuations and the induced charge fluctuations. This gives 
rise to a charge fluctuation operators 



r am (*,E)^n(T,E + hu)at m (E)a 0n (E + ?w). (41) j5*J*( r ,w) = O a/3 (r,w) - /drTI(r,r> Q/3 (r>) 



This equation defines a density of states matrix ^ 7 (r) 
with elements 

v^ Smn {Y,E,E + hw) = h- 1 (v im (E)v Sn (E + hoj)y 1/2 

TP* jm (*,E)^ Sn (r,E + tuj). (42) 

If we form a vector a Q with the annihilation opertors of 
the (incoming) channels in contact a a the local density 
operator in the low frequency limit can be expressed with 
the help of this matrix in the following way, 



n(r, ui) 



a(3 



J dEhi{E)v aP (r)kp{E + wuj). 



(43) 



It is now very convenient and instructive to consider an 
expression for the density operator not in terms of wave 
functions but more directly in terms of the scattering ma- 
trix. In the zero- frequency limit, Ref. ^ introduces the 
matrix 



i//9 7 (a,r) =—t-:[b{ 



4ni l a ^eSU(r) eSU(r) 



(44) 



All quantities in this expression are evaluated at the 
Fermi energy E. The matrix elements of Eq. ( f±2| ) are 
connected to ihe matrices f/3 7 (a, r) (and thus to wave 
functions) viaQ 



E 



7T7 ^ ^ 



2m 



'a/3 



5s 
eSU(r) 



(45) 



Eq. (|45|) is a local Wigner-Smith life-time matrixQ. 

Eq. ( flBI) was given in Ref. |^ and a derivation of this 
relation with the help of an optical potential is given in 
detail in Ref. The relation between scattering states 
and the Green's function of the Schrocdingcr equation on 
one hand and its connection to thi? scattering matrix on 
the othep_hand can also be usedoeJ. Recently Schome- 
rus et al.Ej proposed still another way to proceed by using 
not the overall scattering matrix but scattering matrices 
for different portions of the system. 

It is important to realize that in an electrical conduc- 
tor charge fluctuations are not free but are screened like 



(47) 



Eq. ( |4^ ) permits a gauge invariant dsicussion of charge 
fluctuations in mesoscopic conductors. Eq. (47) has been 
used in Ref. ^ to obtain the second order in frequency 
term of the ac-conductance (see Eq. (|3q)). For a capaci- 
tive structure this defines a novel resistance which deter- 
mines the dissipative effects in charging and decharging a 
mesoscopic conductor. Pedersen and van Langen and the 
authorc3 investigated the current induced into a nearby 
gate due to charge fluctuations in quantum point contacts 
and chaotic cavities. More recently, the charge fluctua- 
tions in two nearby mesoscopic conductorsc3 was treated 
and the effect of quantum dephasing due to charge fluc- 
tuations was calculated within thiS|-aporoachaE3. Often 
self-consistent effects are neglectedcBE]. Charge correla- 
tions at two contacts capacitively coupled to a mesoscopic 
structure are investigated in Ref. |^|. 

Integration of Eq. (45) over the volume of the sam- 
ple gives the global Wigner-Smith life-time matrix. In 
contrast to the local Wigner-Smith matrix the proper- 
ties of the (global) Wigner-Smith life-time matrix have 
received wider attention. In recent years the focus has 
been on the calculation of the entire distribution func- 
tion of delay timescll for structures whose dynamics is in 
the classical limit-chaotic (chaotic cavities) and of one- 
dimensional wiresHoo. While in many cases an analy- 
sis of the global density is quite sufficient, in disordered 
structures, since screening is local, it is clearly desirable 
to proceed from the local density and to investigate the 
properties of the local Wigner-Smith matrix given here. 



XI. DISCUSSION 

We have discussed a hierarchy of density of states as 
they occur in open multiprobe mesoscopic conductors. 
We have illustrated that a small absorption or a small 
emission of particles leads naturally to the appearance of 
these density of states. We have shown that the transmis- 
sion probabilities through weakly coupled contacts, like 
the STM, is related to these densities. We have shown 
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that a weakly coupled voltage probe, describing inelas- 
tic scattering or a dephasing process can be described in 
terms of these densities. We have given a formulation 
of the Friedel sum rule which holds even when the con- 
ditions for the WKB approximation do not apply. We 
have also pointed out that the ac- conductance of a meso- 
scopic conductor at small frequencies can be formulated 
in terms of these densities of states. Furthermore, we 
have indicated that it is useful to consider also the off- 
diagonal elements of a partial density of states matrix 
since this permits a description of fluctuation processes. 
This consideration has led us to a local Wigner-Smith 
life-time matrix. We have emphasized that for many 
problems in electrical conductors, it is necessary to con- 
sider the effect of screening. We have treated screening 
in a simple one-loop random phase approximation. 

The partial density of states which we have used here 
provide a classification with respect to the contact (and 
or channel) index of the scattering matrix. In problems 
with high spatial symmetry^] we can envision also a clas- 
sification with respect to the local momentum. Interest 
in spin transport will require spin resolved partial density 
of states. In hybrid normal superconducting systemstll a 
resolution of both electron and hole density of states is 
conceptually usefulE3e3. 

The description of electrical conduction processes in 
terms of transmission and reflection probabilities has be- 
come very well known. The fact that we can also express 
densities with the help of scattering matrices makes the 
scattering approach applicable not only to dc-transport 
but to a much wider range of electrical transport phe- 
nomena. 
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